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Abstract 

The  concept  being  designed  by  LSRRI  of  united  positioning  and  timing  service  on  the  basis  of  the  uti¬ 
lization  of  long-range  and  global  radionavigation  and  eommon  time  systems  and  aids  for  different  users  is 
deseribed.  The  estimate  of  proposals  of  its  utilization  on  the  national — as  well  as  on  international — scale  is 
given. 


1.  INTRODUCTION 

The  main  direction  of  Leningrad  Scientific  Research  Radio  Technical  Institute  (IjSRRI) 
activities  is  the  designing  and  improvement  of  position  location  and  time  synchronization 
systems  and  aids  for  different  branches  of  national  economy  as  well  as  for  scientific  research. 

The  general  principle  of  systems  and  aids  to  be  developed  by  LSRRI  is  a  generation 
and/or  utilization  of  high-precision  position  and  time  data  for  long-range  and  global  ra¬ 
dionavigation,  .synchronization  of  remote  clocks  and  other  goals  of  this  kind. 

The  most  significant  directions  of  LSRRTs  activities  comprise; 

•  designing  of  equipment  for  transmitting  and  monitor-arid 

•  correction  stations  and  other  aids  of  ground-based  long  range  and  global  pulse-phase 
and  phase  radionavigalion  systems; 

•  creation  of  aids  for  State  System  of  Common  Time  and  Standard  Frequencies,  includ¬ 
ing  monitoring  sites  of  different  accuracy  classes  for  specific  users; 

•  development  of  quartz  oscillators  and  atomic,  frequency  standards  with  optical  (laser) 
pumping  on  the  basis  of  atomic  beam  discriminators  and  hydrogen  masers; 

•  designing  of  on-board  equipment  and  ground-based  aids  for  GLONASS  Satellite  ra¬ 
dionavigation  system  (SRNS); 

•  development  of  position-timing  user  equipment  operating  via  ground-based-as  well 
as  space  based  systems  for  sea,  airspace  and  land  vehicles; 

•  realization  of  scientific  researches. 
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2.  STATE-OF-THE-ART  OF  SOVIET  SYSTEMS  AND 
AIDS  PROVIDING  POSITIONING  AND 

TIMING  SERVICE 

up  to  date,  the  task  of  realization  of  the  positioning  and  timing  service  with  desired 
performances  for  different  platforms  was  being  solved  in  the  USSR,  on  the  basis  of  a  separate 
development  of  aids  for  position  determinations  and  that  for  time  determinations.  The 
those  conditions,  for  position  determinations  was  the  navigation  equipment  being  designed, 
and  for  time  determinations  was  the  timing  equipment  being  developed.  As  much  as  no 
one  of  systems  to  be  mentioned  above  enables  the  achievement  of  desired  positioning  and 
timing  service,  this  task  is  solved  on  the  basis  of  utilization  of  a  system  mix  to  be  optimal 
for  any  region  over  a  certain  time  interval.  As  far  as  we  know,  the  US  specialists  are  of  the 
same  opinion.  It  is  proven  by  US  Federal  Radionavigation  Plan. 

While  developing  the  similar  program  for  long-term  improvement  of  radionavigation 
systems  and  aids  in  the  USSR,  LSRRI  had  already  in  early  ’70  years  has  proposed  idea 
about  a  generation  of  common  global  radionavigation  field  as  a  source  of  the  universal 
position-timing  data  for  any  user  with  the  achievement  of  desired  values  from  the  stand¬ 
point  of  accuracy,  availability,  reliability  and  integrity.  This  approach  was  developed  on 
implementation  of  the  concept  for  the  common  positioning  and  timing  service.  This  con¬ 
cept  has  later  been  approved  and  enlarged  by  a  scries  of  scientific  research  institutes  of  the 
USSR  Academy  of  Sciences  as  well  as  by  ones  of  different  branches  of  national  economy. 

3.  MAIN  CONCEPT  PRINCIPLES  (SUBSTANTIATION 
OF  METHODOLOGY) 

The  basis  of  our  concept  of  providing  a  joint  positioning  and  timing  service  is  a  community 
of  principles  of  position  (coordinates)  and  time  determinations  and,  as  a  consequence,  a 
possibility  to  treat  these  determinations  as  being  two  sides  of  a  single  physical  process  to 
be  connected  with  signal  transfer  and  generation  of  corresponding  data.  As  an  example, 
one  can  use  the  main  navigation  equations  connecting  the  increments  of  space  coordinates 
with  time  increments  to  be  correspondent  with  signal  transfer  between  the  points  to  be 
compared.  In  such  a  case,  navigation  data  are  determined  as  being  data  referring  to  relative 
coordinates  of  bodies  (i.e.  sets  of  material  points)  within  a  four-dimensional  space-time 
variety.  In  its  turn,  the  relative  coordinates  are  the  functions  of  mass  center  mutual  location 
and  mutual  orientation  of  corresponding  axes  a.s  well  as  of  its  derivatives  of  different  degrees. 
Under  these  conditions,  in  order  to  determine  the  space-time  interval  setting  a  relative 
position  of  four  dimensional  points  or,  in  other  words,  to  compare  the  events  to  be  taken 
place  at  those  points,  it  is  necessary  to  use  a  material  data  bearer  moving  within  a  four 
dimensional  framework.  While  providing  the  data  relating  to  origins  of  frames  of  reference, 
the  radio  signals  appear  to  be  such  a  bearer.  The  other  performances  of  navigation-time 
data  depend  on  metrical  properties  of  frame  wherein  the  signal  moves.  Such  an  approach 
enables  to  formulate  a  general  definition  of  any  positioning  and  timing  service  as  a  common 
set  of  space-remote  technical  aids  and  natural  space  sources  to  be  used  for  generation  of 
corresponding  data. 
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4.  FUTURE  NATIONAL-WIDE  POSITIONING  AND 
TIMING  SERVICE  (PRINCIPLES  OF  ORGANIZATION 
AND  OPERATION) 

The  Unified  State  System  of  Positioning  and  Timing  Service  (USS  PTS)  will  be  imple¬ 
mented  on  the  basis  of  existing  and  deploying  ground-and  space-based  navigation  and 
timing  systems  and  aids  by  means  of  institutional  and  methodical  coordination  of  opera¬ 
tion  of  such  aids  with  application  of  equipment  to  be  designed  for  solving  the  user  problems. 
This  is  equally  valid  for  the  equipment  of  Joint  position-  timing  complexes. 

A  specific  feature  of  such  a  system  is  its  openness,  permitting  a  possibility  of  including 
as  well  as  exclusion  of  certain  aids.  The  other  important  feature  is  a  possibility  of  generation 
of  unified  positioning-timing  basis  providing  a  performing  of  determinations  in  the  unified 
frame  of  reference  with  a  possible  gradation  in  accordance  with  accuracy  classes  and  with 
a  subsequent  transition  to  any  frame  of  reference  to  be  convenient. 

From  the  organizational  standpoint,  USS  PTS  might  consist  of  four  subsystems  as 
follows; 

•  fundamental  bases  of  space  and  time  frame  of  reference; 

•  generation  and  transfer  of  position-timing  data; 

•  control  and  monitoring  of  system  operation; 

•  users. 

The  important  feature  of  the  system  to  be  proposed  is  that  it  isn’t  merely  a  set  of 
technical  aids  or  an  organization  operating  in  accordance  with  the  set  of  rules  to  be  pre¬ 
scribed,  but  that  is  a  techidcal  system  of  a  top  level.  Such  a  system  may  be  considered  as 
a  ^synergetic  one  permitting  to  realize  a  new  quality-the  presentation  of  space-time  data 
with  ensuring  the  values  of  accuracy,  availability,  reliability  and  integrity  to  be  necessary 
for  defined  users  while  solving  all  the  problems,  including  specific  ones. 

The  USS  PTS  implementation  on  the  basis  of  RNS  to  be  already  in  operation  and 
to  be  deployed  would  enable  to  set  in  foreseeable  future  a  global  positioning-time  field 
with  accuracy  level  of  10m  (Im  locally)  and  synchronization  error  of  50ns,  thus  providing 
a  possibility  of  solving  the  widest  class  of  problems  for  all  users,  including  a  short-range 
navigation  task  with  landing  approach  and  landing  itself. 

5.  PROSPECTS  OF  INTERNATIONAL  COOPERATION 
WITHIN  THE  FIELD  OF  GLOBAL  PTS  CREATION 

From  our  point  of  view,  the  national-wide  principles  of  PTS  organization  to  be  considered 

above  would  be  used  as  the  basis  of  concept  of  global  PTS  implementation  - .  of  course, 

with  taking  into  account  a  number  of  factors  relating  to  foreign  political  situation,  USSR 
and  USA  national  priority  systems,  specific  roles  of  certain  regions  as  well  as  existing 
machineries  for  attainment  of  agreements.  The  advantages  of  an  open  positioning  and 
timing  service  may  be  realized  only  in  conditions  of  “open”  international  community,  i.c. 
such  a  community  wherein  the  relations  between  the  countries-participants  are  based  on 
the  whole  mutual  confidence.  With  regard  to  this,  how  can  one  estimate  the  relations 
between  the  superpowers,  the  USSR  and  the  USA,  and  its  influence  on  prospects  of  global 
PTS  creation? 
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With  the  beginning  of  current  decade,  the  world  to  come  in  the  initial  period  of  the 
new  non-confrontation  era,  but  this  period  is  complicated  and  contradictory,  because  the 
past  tries  to  keep  its  positions  and  the  future  carries  the  non-predictable  dangers.  The 
nuclear  arsenals  lose  its  sense,  but  its  improvement  is  going  on.  The  space-based  PTS 
aid  that  appear  to  be  intended  for  the  main  rriole  on  future  PTS  implementation  on  a 
national-wide  as  well  as  on  international-wide  scale  were  being  created  in  the  USSR  and 
in  the  USA,  in  conditions  of  competition  for  the  defense  agencies  with  taking  into  account 
its  tasks  in  the  first  place.  1’he  very  concept  of  defense  sufficiency  gives  today  only  certain 
quantity  cut  of  armaments,  while  not  providing  an  ending  of  its  quality  improvement.  It 
is  equally  true  for  the  supporting  subsystems,  including  PTS  aids. 

In  such  a  situation,  the  only  possible  way  to  overcome  a  distrust  to  be  accumulated 
over  the  preceding  decades  of  the  “Cold  War”  is  that  to  be  based  on  an  improvement  of  a 
machinery  for  generation  and  adoption  of  .Toint  solutions. 

From  our  point  of  view,  it  should  prefer  any  form  of  cooperation,  including  ramified 
membership  in  different  institutions,  creation  of  cross  groups  and  parallel  institutes,  partly 
duplicating  its  functions,  but  permitting  at  the  same  time  to  achieve  the  best  mutual 
understanding  as  a  result. 

At  the  present  time,  LSRRI’s  practical  activities  include  examples  of  cooperation  with 
the  USA  and  Western  European  state  agencies  as  well  as  with  different  companies. 

The  development  and  the  substantiation  of  a  concept  with  taking  into  account  the  de¬ 
mands  of  different  branches  of  national  economy  as  well  as  the  needs  of  scientific  researches 
are  being  carried  on  LS  RRl  over  a  series  of  years.  In  November  of  1990,  in  Leningrad,  the 
Fist  All-Union  Scientific  and  Technical  Conference  on  problems  relating  to  the  common 
positioning  and  timing  service  was  held  at  LSRRL  It  is  planned  for  such  conferences  to  be 
the  regular  ones  in  future. 

In  conclusion,  I  would  like  to  note  that  there  is  a  sufficiently  powerful  factor  permitting 
an  optimistic  estimate  for  the  prospects  of  realization  of  global  PTS.  This  factor  is  a 
belonging  of  such  a  system  realizing  the  needs  of  the  civilized  community  on  a  global 
scale  to  the  values  common  to  humanity  and,  consequently,  its  incontestable  priority.  The 
support  of  a  safe  operation  of  such  a  system  will  undoubtedly  need  the  highest  art  not 
only  with  regard  to  scientific  technical  sphere,  but  also  to  human  contacts.  We  would  like 
to  hope  the  good  will  of  all  participants  would  enable  a  successful  solving  of  this  problem. 
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GPS  COMPOSITE  CLOCK  SOFTWARE  PERFORMANCE 

A.  L.  Satin,  W  A.  Feess,  H.  F.  Fliegel,  and  C.  H.  Yinger 
The  Aerospace  Corporation 
Los  Angeles,  CA  90009 

ABSTRACT 

Computer  software  for  ensembling  all  the  space  vehicle  (SV)  and  ground  clocks  of  the  Global  Positioning 
Systems  (GPS)  was  implemented  at  the  Master  Control  Station  (MCS)  on  17  June  1990.  Improved  GPS 
time  stability  and  steering  control  are  predicted.  This  paper  assesses  the  perceptions  of  both  Air  Force  and 
outside  users  and  compares  current  performance  to  theoretical  predictions. 

GPS  SYSTEM  DESCRIPTION 

The  Global  Positioning  System  is  a  Department  of  Defense  space-based  navigation  and  timing  dis¬ 
semination  system  currently  under  deployment.  By  1993  it  will  consist  of  a  constellation  of  24  space 
vehicles  (SVs)  supported  by  a  control  segment  (CS)  that  tracks  the  SVs  with  its  monitor  stations 
(MSs),  determines  SV  ephemeris  states  and  SV  and  MS  clock  states  at  the  Master  Control  Station 
(MCS)  and  periodically  uploads  the  SVs  with  predicted  SV  states  that  users  require  to  navigate  and/ 
or  determine  their  time.  At  this  stage  in  the  buildup,  the  constellation  consists  of  six  block  I  SVs  and 
nine  block  II  SVs.  The  CS  MSs  are  located  at  Kwajalein  Island,  Diego  Garcia  Island,  Hawaii,  Ascen¬ 
sion  Island,  and  Colorado  Springs,  Colorado. 

SV  range  measurements  (i.e.,  pseudorange)  are  a  function  of  the  SV(i)-MS(j)  geometric  range  and  the 
phase  offset  between  the  SV(i)  and  MS(j)  clocks,  where  i  and  j  are  the  ith  SV  and  the  jth  MS,  respec¬ 
tively.  The  accuracy  with  which  a  user  can  navigate  (i.e.,  determine  his  position  and  clock  state)  de¬ 
pends  on  the  accuracy  of  his  measured  pseudorange  and  the  accuracy  of  the  predicted  SV  ephemeris 
and  clock  information,  generated  and  uploaded  by  the  CS.  The  latter  source  of  error  is  fundamental 
to  GPS  operation  and  is  quantified  by  the  User  Range  Error  (URE).  URE  is  defined  as  the  differ¬ 
ence  between  the  range-to-a-user  computed  using  predicted  states  and  the  “true”  range-to-a-user 
computed  using  the  latest  Kalman  estimates  of  ephemeris  and  clock.  Operationally,  the  MCS  com¬ 
putes  UREs  every  15  min  (the  Kalman  filter  update  cycle)  for  each  SV,  taking  a  root  sum  of  squares 
(RMS)  of  individual  UREs  to  a  fixed  set  of  grid  locations  on  the  earth  that  have  visibility  to  the  SV. 

The  CS,  by  solving  for  SV  and  MS  clock  phase  and  frequency  states  from  the  pseudorange  measure¬ 
ments,  effectively  synchronizes  the  system  to  GPS  time.  Originally  the  GPS  time  reference  was  taken 
to  be  one  of  the  MS  clocks.  States  for  that  MS  (i.e.,  the  master  MS)  were  defined  to  be  fixed  at  zero 
and,  hence,  were  not  included  in  the  filter  state  vector.  Clock  state  offsets  estimated  by  the  filter  were 
then  relative  to  master  clock  (GPS)  time.  When  the  filter  is  allowed  to  estimate  all  clock  states  (i.e., 
including  states  for  what  used  to  be  the  master  MS),  the  configuration  is  called  a  composite  or  en¬ 
semble  clock  system.  Since  June  of  this  year,  the  time  reference  for  the  system  clocks  (i.e.,  GPS  time) 
has  been  a  composite  clock. 
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Since  the  MCS  has  so  many  SV  states  to  estimate  (11  per  SV  x  24)  it  was  decided  long  ago  to  “parti¬ 
tion”  the  SVs  among  several  filters  each  capable  of  handling  up  to  six  of  them  and  each  containing,  of 
course,  the  complete  complement  of  MS  states:  3  per  MS  x  4  for  a  master  clock  configuration  and  3 
per  MS  X  5  for  the  composite  clock  case. 

COMPOSITE  CLOCK  IMPLEMENTATION  AND  THEORY 
When  additional  states  are  added  to  the  Kalman  filter,  for  what  was  the  master  MS,  the  GPS  system 
becomes  explicitly  unobservable.  This  is  because  a  constant  bias  shift  in  all  clock  phase  states  would 
not  affect  pseudorange  and  hence,  would  be  unobservable  to  the  filter.  The  existence  of  this  unob¬ 
servable  component  of  state  causes  the  clock  phase  covariance  matrix  to  grow  linearly  with  time  even 
though  the  system  is  completely  stable  with  well  behaved  filter  gains.  See  Reference  1  for  further  dis¬ 
cussion  of  these  characteristics.  Since  a  growing  covariance  matrix  will  eventually  cause  numerical 
problems,  a  method  for  reducing  the  covariance  matrix  without  degrading  filter  performance  was  in¬ 
troduced.  The  method  chosen  was  reduction  via  pseudomeasurement  update  (Reference  2). 

Covariance  Reduction 

This  update  takes  the  standard  Kalman  form 

F’  =  P-lPtfiHFff  +  R)-']  HP  (1) 

where  P'  is  the  covariance  matrix  resulting  from  the  update  (i.e.,  the  reduced  covariance  matrix)  and 
R  is  the  pseudomeasurement  noise  variance.  The  matrix  P  in  Equation  1  is  the  full  n  x  n  covariance 
matrix  of  ephemeris  and  clock  states  for  the  SVs  and  MSs  of  a  particular  partition.  All  filter  parti¬ 
tions  undergo  independent  covariance  reductions  every  filter  update  cycle.  The  n  x  1  column  vector 

is  constructed  by  inserting  each  element  of  an  m  x  1  column  vector  into  the  appropriate  posi¬ 
tion  in  an  n  x  1  vector  of  zeros.  is  given  by 

(2) 

where  B  is  the  m  X  m  submatrix  of  P  corresponding  to  clock  phase  states  only  and  u  is  an  m  x  1  col¬ 
umn  vector  of  elements,  all  of  which  are  one.  Each  clock  in  a  partition,  then,  has  an  ensembling 
weight  in  the  column  vector  H  (ordered  as  P  is  ordered),  whereas  each  ephemeris  state  has  weight 
zero.  By  the  construction  of  H,  the  sum  of  all  weights  is  one.  Equation  2  assigns  to  each  clock  a  nor¬ 
malized  set  of  phase  and  frequency  weights  that  are  inversely  proportional  to  the  magnitude  of  esti¬ 
mation  error  variance  (i.e.,  the  diagonals  of  B).  Because  of  the  particular  form  of  H^,  however,  the 
correction  term  in  Equation  1,  i.e.,  the  second  term  on  the  right  side  of  the  equation  has  nonzero  ele¬ 
ments  corresponding  to  ephemeris  states  and,  hence,  this  method  of  covariance  reduction  is  not  a 
transparent  variation  on  filter  operation.  A  simulation  of  this  algorithm  was  performed  (Reference  3) 
using  real  pseudorange  measurement  data  from  the  MCS  in  order  to  assess  the  impact  of  the  algo¬ 
rithm  on  Kalman  states  and  URE.  RMS  UREs  were  shown  to  be  very  good  with  the  algorithm  and 
states  changed  only  slightly.  Ephemeris  states  with  composite  clock  operation  are  unchanged  from 
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what  they  were  with  the  master  clock  configuration— with  or  without  the  reduction  algorithm  (Refer¬ 
ence  3). 

Partition  Reconciliation 

Partition  reconciliation  is  required  to  ensure  that  composite  clock  time  for  the  various  filter  partitions 
is  consistent.  Uploads  built  for  SVs  from  different  partitions  must  have  very  nearly  the  same  refer¬ 
ence  time  for  clock  phase  and  frequency  offsets  contained  in  the  navigation  message.  Any  difference 
in  reference  times  maps  into  URE.  The  reconciliation  process,  since  it  is  an  adjustment  of  reference 
time  for  the  partition,  does  not  affect  filter  gains  or  measurement  residuals  because  it  amounts  to  a 
constant  shift  of  all  the  clock  phase  and  frequency  states  in  the  partition.  Each  partition  will  have  a 
different  shift  as  computed  by  the  reconciliation  algorithm.  Partition  reconciliation  is  possible  only 
because  the  MS  clocks  are  common  to  each  partition.  It  requires  that  a  specified  weighted  sum  of 
MS  states  be  the  same  for  each  partition,  where  the  same  weights  are  used  to  compute  the  sum  for 
each  partition.  The  value  of  the  weighted  sum  is  itself  a  weighted  average  of  MS  and  SV  states.  The 
computation  proceeds  as  follows: 

First  the  vectors  b,  are  computed  for  each  partition  according  to  (References  2,  4) 


b,  =  [(//  *]-\H  *fC-^ 


(3) 


where  ^  are  the  MSj  clock  states  for  the  i*  partition  and  H*  is  a  10  x  2  matrix  of  five  2x2  identity 
matrices.  C  is  a  10  x  10  matrix  obtained  by  summing  the  MS  clock  covariance  submatrices  from  the 
various  partitions.  Note  from  Equation  3  that  the  linear  combination  of  MS  states  is  taken  for  each 
partition.  Next,  the  vector  b  that  the  bj  are  reconciled  to,  is  computed  according  to  Equation  4  below: 

i“l 


where 


N  =  the  number  of  partitions,  and  nj  =  the  number  of  SVs  in  partition  i. 


j  ^  +  9ir^/3  <7ir^/2 

y  2  /T  J 

<72^72  ^/2r 

for  the  J**’  SV  in  partition  i.  In  Equation  5,  2/,  and  efi  arc  the  phase  and  frequency  process  noise  co- 
variances,  respectively,  and  t  is  the  update  interval  (15  min).  The  adjustment  (Abj)  that  gets  applied 
to  each  clock  phase  state  in  partition  i  is  then  computed  by 


A*,  =^b-bi 


531 


After  the  adjustments  are  made  to  each  partition,  a  recomputation  of  the  bj  would  yield  b  for  all  par¬ 
titions.  This  means  that  GPS  time  for  each  partition  is  now  the  same,  i.e.,  after  reconciliation  the 
weighted  average  of  MS  states  for  each  partition  is  the  same. 

Stability  of  GPS  Time 

One  of  the  primary  reasons  for  ensembling  the  GPS  clocks  (i.e.,  constructing  a  composite  clock)  is  to 
gain  stability  in  GPS  time.  Clock  ensembling  is  a  proven  technology  most  prominently  used  by  the 
U.S.  Naval  Observatory  (USNO)  and  National  Institute  of  Standards  and  Technology  (NIST).  The 

enhancement  of  stability  for  a  system  of  n  clocks,  each  of  which  is  stable  to  say  IE- 13,  goes  as  l/Jii . 
For  the  current  system  of  19  clocks,  operating  with  a  Kalman  filter  configured  for  the  composite  clock 
option,  the  stability  of  GPS  time  is  predicted  to  be  of  the  order  of  2.3E-14.  This  assumes  that  the 
three  MCS  filter  partitions  (with  5,  5  and  4  SVs  in  partitions  1,  2,  and  3,  respectively)  are  being  recon¬ 
ciled  perfectly  and  states  do  not  differ  appreciably  from  what  they  would  be  if  the  MCS  were  operat¬ 
ing  with  the  14  SVs  in  one  partition. 

Steering  of  GPS  Time 

The  composite  clock  software  package  provides  a  new  GPS  time  steering  algorithm  that  automatically 
determines  and  applies  the  appropriate  clock  drift  rate  steer  to  all  clock  states  when  enabled  by  oper¬ 
ator  directive.  The  operator  need  only  input  daily  USNO  updates  of  GPS-UTC  phase  offset.  A 
bang-bang  controller  (Reference  4)  has  been  implemented.  Allowable  steering  rates  are  ±  2E-19s/s^ 
and  0.  Three  steering  modes  are  provided  as  follows:  steer  out  the  GPS-UTC  offset  as  computed  by 
the  MCS  from  “raw”  USNO  measurements;  steer  out  a  particular  GPS-UTC  offset  input  by  the  oper¬ 
ator;  or  steer  GPS  time  to  one  of  the  MS  clocks.  Simulations  of  this  control  loop,  carried  out  by  The 
Aerospace  Corporation  (Reference  5),  predict  an  RMS  steady  state  control  error  of  the  order  of  10  ns 
assuming  the  GPS  time  stability  is  3E--14s/s  at  1  day  (Figure  1). 

COMPOSITE  CLOCK  PERFORMANCE  (ACTUAL) 

Transition  to  the  composite  clock  was  initiated  on  17  June  1990  by  including  the  master  MS  clock 
States  as  estimated  states  in  the  filter  partitions.  This  act  causes  the  state  component  covariances  to 
then  be  Q’d  (i.e.,  process  noise  added  in  the  filter  propagation  operation)  every  15  min. 

UREs 

RMS  UREs  for  all  age-of-data  were  as  good,  if  not  better,  than  what  they  were  with  the  master  clock 
configuration.  For  one  upload  per  day  per  SV  these  RMS  values  were  in  the  range  of  2  to  4  m.  Note 
that  UREs  are  very  insensitive  to  GPS  time  stability.  They  are  primarily  due  to  SV  clock  phase/fre¬ 
quency  variations  over  the  upload  interval. 


At  23:30  on  25  June  1990  (MJD48067),  a  directive  was  input  to  the  system  to  begin  steering  GPS  time 
to  drive  the  GPS-UTC  internally  computed  offset  to  zero.  Figure  2  shows  a  history  of  GPS-UTC  ob¬ 
tained  from  the  USNO  GPS VI  file.  The  steering  achieved  is  not  as  good  as  clock  ensembling  can 
theoretically  achieve.  Control  is  not  able  to  reduce  the  large  amplitude  ( ±  100  ns)  of  the  offset.  It  is 
apparent  that  the  underlying  stability  of  GPS  time  (without  steering)  must  be  considerably  worse  than 
2.3E-14  and/or  some  other  unknown  type  of  aging-like  effect,  perhaps  caused  by  system  model  errors, 
must  be  going  on.  Table  1  shows  the  time  steer  history  from  0000  6  January  1980  to  0630  22  October 
1990.  Note  that  steer  control  values  (i.e.,  ±2.00E-19  or  O.OOE-l-00)  only  change  at  times  that  are  mul¬ 
tiples  of  15  min.  (the  Kalman  filter  update  interval). 

Table  1.  Time  Steer  History 


Day 

Month 

Year 

Modified 
Julian  Date 

Time 

Steer 

Control  Value 

6 

Jan 

1980 

44244 

00:00:00 

+  0.00E  +  00 

25 

Jun 

1990 

48067 

23:30:00 

+  2.00E-19 

11 

Jul 

1990 

48083 

22:30:00 

-2.00E-19 

12 

Jul 

1990 

48084 

15:00:00 

+  2.00E-19 

13 

Jul 

1990 

48085 

15:00:00 

-2.00E-19 

13 

Jul 

1990 

48085 

21:45:00 

+  2.00E-19 

16 

Jul 

1990 

48088 

01:15:00 

-2.00E-19 

16 

Jul 

1990 

48088 

18:30:00 

+  2.00E-19 

17 

Jul 

1990 

48089 

05:15:00 

-2.00E-19 

17 

Jul 

1990 

48089 

17:30:00 

-I-2.00E-19 

19 

Jul 

1990 

48091 

17:00:00 

-2.00E-19 

17 

Aug 

1990 

48120 

16:15:00 

-l-O.OOE-l-00 

27 

Aug 

1990 

48130 

23:45:00 

-2.00E-19 

28 

Aug 

1990 

48131 

06:30:00 

+  2.00E-19 

28 

Aug 

1990 

48131 

17:00:00 

-2.00E-19 

29 

Aug 

1990 

48132 

09:30:00 

+  2.00E-19 

29 

Aug 

1990 

48132 

22:15:00 

-l-O.OOE  +  00 

4 

Sep 

mmn 

48138 

16:45:00 

+  2.00E-19 

17 

Sep 

BBI 

48151 

14:15:00 

+  0.00E  +  00 

21 

Sep 

1990 

48155 

12:15:00 

+  2.00E-19 

21 

Sep 

1990 

48155 

12:45:00 

-2.00E-19 

5 

Oct 

1990 

48169 

18:30:00 

+  0.00E-I-00 

9 

Oct 

1990 

48173 

07:15:00 

-2.00E-19 

Oct 

1990 

48184 

11:15:00 

+  0.00E-I-00 

Oct  1990  48186  06:30:00 


When  the  steering  from  Table  1  is  backed  out  of  the  GPS-UTC  plot  of  Figure  2,  along  with  a  mean 
slope  of  4  ns/day,  the  underlying  behavior  of  GPS-UTC  is  evident  (Figure  3).  Ideally  Figure  3  should 
look  like  random  walk  with  a  small  amplitude  starting  from  250  ns  offset.  The  oscillatory  nature  of 
this  curve  needs  to  be  understood.  The  only  source  of  perturbation  on  GPS  time  (or  GPS-UTC)  must 
come  via  the  Kalman  filter  and  will  occur  when  the  filter  is  “lied  to.”  For  example,  if  a  Rubidium 
clock  (e.g.,  SV  16)  has  a  large  frequency  jump  and  the  clock  frequency  state  is  not  Q’d  sufficiently  to 
absorb  the  variation,  then  GPS  time  will  be  perturbed.  It  would  also  be  perturbed  if  an  MS  suffers  a 
phase  or  frequency  step  that  goes  unnoticed  and  so  the  clock  covariance  is  never  Q-bumped  to  allow 
the  filter  to  solve  for  the  new  phase  or  frequency.  GPS  time  is  perturbed  because  the  corrected  clock 
(i.e.,  the  clock  time  corrected  by  the  filter  state  offset)  represents  GPS  time,  and  if  the  state  doesn’t 
correctly  compensate  for  true  clock  behavior,  the  “corrected”  clock,  or  GPS  time  will  be  perturbed. 

In  the  course  of  these  investigations  perturbations  in  both  MS  and  SV  clock  states  were  evident,  and 
it  is  only  necessary  to  see  how  the  filter  was  handled  (e.g.,  was  the  covariance  matrix  Q-bumped  or 
process  noise  increased  at  these  times)  in  order  to  assess  whether  or  not  a  perturbation  in  GPS  time 
was  likely.  An  attempt  was  made  to  correlate  state  perturbations  occurring  in  the  interval  2  October 
to  15  October  with  the  slope  discontinuities  of  Figure  3.  On  2  October  the  Multiple  Frequency  Stan¬ 
dard  System  (i.e.,  the  MFSS  of  six  Cesium  clocks  operated  by  the  U.  S.  Naval  Research  Laboratory 
(NRL))  at  Falcon  Air  Force  Base  became  the  primary  standard  for  the  Colorado  Springs  Monitor 
Station  (COSPM).  At  the  time  of  the  switch,  a  large  jump  in  the  COSPM  frequency  state  was  appar¬ 
ent.  Hopefully  the  jump  was  simultaneous  with  a  large  Q-bump  of  the  COSPM  clock  states,  however, 
the  distinct  change  in  slope  of  GPS-UTC  at  that  time  is  evidence  to  the  contrary. 

SV  Clock  Stability 

Although  GPS  time  is  being  perturbed,  due  to  unmodeled  filter  processing,  the  stability  of  the  SV 
clocks  is  up  to  expectations.  Figure  4,  provided  by  D.  Allan  of  NIST  is  a  plot  of  Allan’s  variance  for 
five  SVs  in  the  constellation.  SV  16,  the  outlier,  was  having  problems  during  this  interval  of  time  (12 
August  to  11  September  1990).  The  other  four  block  II  SVs  are  operating  better  than  their  specifica¬ 
tion  of  2E-13  at  tau  equals  1  day.  Figures  5  and  6  show  SVN  18  clock  phase  behavior  (with  a  straight 
line  removed)  and  corresponding  Allan  variance  for  tau  between  15  min  and  1  day.  These  were  com¬ 
puted  from  MCS  Kalman  filter  data.  A  12-hr,  unmodeled,  periodic  variation  in  phase  is  evident  in 
Figure  5.  It’s  peak-to-peak  amplitude  is  about  2  m.  This  periodicity  produces  the  sharp  breaks  in 
the  Allan  variance  plot  of  Figure  6.  The  1-day  value  of  the  variance,  8E-14,  is  uncorrupted  by  the 
periodicity,  since  1  day  is  a  multiple  of  12  hr. 

Stability  of  GPS  Time 

It  should  come  as  no  surprise,  because  of  the  oscillatory  behavior  of  GPS-UTC  evidenced  in  Figure  3, 
that  GPS  stability  relative  to  UTC  (in  the  interval  31  May  to  8  September  1990)  would  not  be  as  good 
as  predicted  from  the  theory  (2.3E-14),  This  is  apparent  from  Figure  7,  another  of  D.  Allan’s  contri¬ 
butions.  Stability  is  more  like  8E-14;  comparable  to  the  performance  of  a  master  MS  configuration. 
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It  should  be  noted  that  in  the  time  span  of  Figure  7,  SV  16  had  a  large  frequency  jump  and  in  addi¬ 
tion  its  aging  state  was  not  being  Q’d  in  the  filter.  This  could  have  had  a  significant  impact  on  GPS 
time.  Another  impression,  however,  of  GPS  time  stability  is  gained  when  one  looks  at  the  Colorado 
Springs  monitor  station  clock  phase  state.  This  state  is  an  estimate  of  the  phase  offset  between  the 
MFSS  ensemble  (stable  to  2E-14  relative  to  UTC)  and  GPS  time.  The  phase  state  from  partition  No. 
1  is  plotted  in  Figure  8  and  its  Allan  variance  is  shown  in  Figure  9.  The  1-day  GPS  stability  (with  the 
MFSS  stability  removed)  is  4E-14.  The  difference  between  this  assessment  (based  on  the  current 
time  phase  state)  and  that  of  D.  Allan  (based  on  the  predicted  phase  state)  is  considerable.  The 
1-day  Allan  variance  is  close  to  our  theoretical  expectation.  Figure  10,  an  Allan  variance  plot  for  the 
Kwajalein  MS,  also  supports  the  idea  that  GPS  time  stability  is  better  than  4E-14  at  1  day.  The  24-hr 
periodicity  inherent  in  the  clock  phase  of  Kwajalein  (and  all  other  MS  clocks  to  a  greater  or  lesser 
degree)  causes  the  Allan  variance  to  be  larger  at  12  and  36  hr  than  at  24  hr.  Figure  11  shows  this  pe¬ 
riodicity  with  peak-to-peak  amplitude  of  4  m  or  so.  If  the  clocks  are  really  cycling  daily  (e.g.,  because 
of  temperature  and/or  humidity  changes),  and  the  filter  estimates  only  reflect  a  portion  of  the  ampli¬ 
tude,  it  is  possible  that  GPS  time  stability  will  be  degraded  for  taus  that  are  not  a  multiple  of  24  hr 
(as  Figures  9  and  10  indicate).  On  the  other  hand,  augmenting  the  clock  model  with  a  sinusoid  term 
would  remove  the  performance  degradation.  It  is  also  possible  that  the  unmodeled  sinusoid  is  hin¬ 
dering  our  ability  to  steer  GPS  time.  This  is  being  investigated  by  simulation. 

As  mentioned  earlier,  each  filter  partition  has  its  own  GPS  time  reference,  and  these  are  reconciled 
every  15  min  for  a  common  value.  The  effectiveness  of  this  reconciliation  is  measured  by  differencing 
corresponding  MS  phase  states  from  two  partitions.  One  of  these  differences  for  the  Kwajalein  MS 
and  partitions  No.l  and  3  is  plotted  in  Figure  12.  Excursions  are  generally  between  ±2m.  The  Allan 
variance  of  the  difference,  a  direct  measure  of  the  stability  of  partition  No.l  GPS  time  relative  to 
partition  No.  3,  is  4E-14  at  1  day  (Figure  13).  The  inferred  absolute  stability  of  each  relative  to  a 
“perfect”  clock  would  be  4/v^  E- 14,  or  2.7E-14.  The  stability  of  GPS  time  for  each  partition  looks 
good  in  terms  of  Allan  variance,  but  the  MS  phase  difference  of  ±  2  m  may  be  a  problem.  This  is 
being  investigated  further. 


CONCLUSIONS 

GPS  time,  constructed  from  the  composite  clock  software,  is  being  perturbed  by  unmodeled  or  im¬ 
properly  modeled  clock  phase  behavior.  Analysis  of  Kalman  filter  MS  clock  phase  estimates  shows  a 
daily  periodic  signature  with  peak-to-peak  amplitude  of  4  m  that  varies  somewhat  from  MS  to  MS. 
GPS  time  stability,  inferred  from  Allan  variances  of  MS  clock  phase  estimates,  in  the  interval  2  to  22 
October  1990  very  nearly  measures  up  to  expectations.  This  means  that  the  filter  clock  states  are  by 
and  large  reflecting  true  clock  behavior  and,  hence,  are  isolating  GPS  time  from  clock  perturbations 
in  this  interval.  Further  studies  are  needed  to  explain  the  behavior  of  GPS-UTC  in  the  interval  prior 
to  2  October  when  SV  16  was  having  its  problems. 
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NANOSECONDS 


47700  47780  47860  47940  48020  48100  48180  48260 


JULIAN  DATE  (note:  MJD  47700  =  23  June  1989) 

Figure  2.  UTC(USNO)  Reference  —  GPS  Time 
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NANOSECONDS 


29  JUNE  THROUGH  29  OCTOBER  1990 


Figure  3.  GPS-UTC(USNO)  with  GPS  Time  Steering  and  Mean  Slope  Removed  (4  ns/day). 
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Figure  4.  PRN(SVN)  Clock  —  UTC(NIST)  12  August  to  11  September  1990. 
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METERS 


TIME  (day  of  October) 

Figure  8.  Colorado  Springs  Monitor  Station  Clock  Phase  Estimate  (Partition  No.  1). 
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METERS  METERS 


TIME  (day  of  October) 


Figure  11.  Kwajalein  Monitor  Station  Clwk  Phase  Estimate  (3-19  October  1990). 


TIME  (day  of  October) 

Figure  12.  Partition  1  Minus  Partition  3  Kwajalein  Monitor  Station  Clock  Phase  Difference. 
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SECONDS 


Figure  13.  Allan  Variance  of  Partition  1  Minus  Partition  3  Kwajalcin  Monitor  Station  Clock  Phase  Difference. 
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QUESTIONS  AND  ANSWERS 


Professor  Alley,  University  of  Maryland:  You  mentioned  that  there  was  a  variation.  Did 
that  have  a  12  hour  or  a  24  hour  period? 

Mr.  Satin:  That  did  not  have  any  definite  periodicity  to  it.  It  does  look  like  the  station  frequency 
has  a  24  hour  periodicity  to  it.  We  think  that  the  periodicity  is  caused  by  humidity. 

Professor  Alley:  You  mentioned  a  12  hour  period  also. 

Mr.  Satin:  That  is  in  the  satellite,  they  have  a  12  hour  orbit.  That  is  a  temperature  variation 
effect. 

Dr.  Winkler,  U.  S.  Naval  Observatory:  What  was  the  rationale  which  led  to  the  decision  to 
use  a  bang-bang  steering  to  UTC?  And  what  is  the  dead-hand  in  that  control? 

Mr.  Satin:  I  can’t  answer  the  question  about  the  dead-band  right  now.  We  investigated  several 
different  algorithms  for  the  control.  The  bang-bang  takes  out  large  initial  offsets  much  quicker  than 
the  others. 

Dr.  Winkler:  At  the  expense  of  much  larger  frequency  changes,  and  that  is,  of  course,  detrimental 
to  its  time  usage. 

Mr.  Satin:  All  the  simulations  show  that  there  is  no  problem  there. 

David  Allan,  NIST:  On  your  last  view-graph,  the  cesium  performance  of  from  10“^^  down 
to  10"^^  is  pretty  typical  of  the  space  cesiums.  How  can  you  have  a  type  of  behavior  out  at  10® 
seconds?  That  is  not  typical  of  cesiums  or  anything  in  the  system. 

Mr.  Satin:  We  think  that  that  is  related  to  the  periodicity  in  the  monitor  stations.  In  other  words, 
at  the  24  hour  point  you  don’t  see  the  24  periodicity  and  things  get  better,  but  you  do  see  it  at  other 
points. 

Mr.  Allan;  You  are  saying  that  this  is  an  aliasing? 

Mr.  Satin:  Yes,  that  is  right. 

Mr.  Allan;  But  if  you  have  a  periodicity,  it  will  increase  the  noise  over  it 

The  curve  would  have  an  increase  at  12  hours,  and  no  increase  at  24  hours.  You 
so  I  don’t  see  any  aliasing  affect. 

Dr.  Henry  Fliegel,  The  Aerospace  Corporation:  We  believe  that  the 
of  the  sigma-tau  curve  is  because  we  do  have  24  hour  periodicities  and  they  alias 
As  far  as  I  know,  that  completely  explains  the  appearance  of  the  curve. 

Dr.  Claudine  Thomas,  BIPM:  Why  are  you  using  that  set  of  filters  and  the  clocks  are  not 


will  not  decrease  it. 
have  a  curve. 


strange  appearance 
at  12  and  36  hours. 
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independent,  when  with  one  Kalman  filter  you  could  handle  everything? 

Mr.  Satin:  It  is  compntationally  more  efficient  to  use  three  smaller  filters.  When  the  system  was 
first  designed,  we  didn’t  have  the  number-crunching  capability  to  handle  the  single  Kalman  filter. 
That  would  be  24  satellites,  11  states  per  satellite  for  a  total  of  hundreds  of  states. 

Dr.  Thomas:  Yes,  but  the  clocks  which  are  involved  in  your  filters  are  not  independent. 

Mr.  Satin:  Right,  it  is  sub-optima,!.  That  is  correct,  but  we  did  a  study  to  determine  what  the 
effect  would  be  to  use  one,  two,  or  three  filters  because  we  are  making  a  large  software  change  and 
could  change  the  filters,  but  the  gain  was  essentially  zero. 

Sam  Stein,  Ball  Corp.:  I  am  afraid,  Henry,  that  I  would  have  to  agree  with  Dave  Allan  and 
disagree  with  you.  We  calculated  the  exact  effect  of  a  periodicity  on  Oy.  Its  effect,  if  you  have  a  24 
hour  periodicity,  would  be  an  elevation  of  ay  over  the  white  frequency  noise  level,  due  to  the  ensemble, 
at  12  hours,  a  decrease  back  to  the  white  frequency  noise  level  at  24  hours,  and  then  it  would  come 
back  up  afterwards.  If  you  look  at  the  very  short  term,  if  you  look  at  the  numbers  at  1000  seconds, 
you  can  see  where  the  level  is.  Whatever  this  phenomenon  is,  it  causes  a  decrease  below  the 

noise  level  that  one  would  attribute  to  the  ensemble,  and  that  is  surprising.  We  often  see  this  when 
there  is  a  small  amount  of  data,  but  if  those  are  error  bars  that  I  see  on  the  curve,  it  would  indicate 
that  there  is  an  adequate  amount  of  data  and  one  should  not  see  such  a  large  effect. 

Mr.  Satin:  We  have  approximately  20  days  of  data,  so  we  should  have  a  pretty  good  estimate  of 
the  one  day  variance. 

Tom  McCaskill,  Naval  Research  Laboratory:  I  would  just  like  to  make  some  comments.  If 
you  look  at  the  short  term  data,  you  start  out  at  900  seconds,  which  is  15  minutes.  If  you  follow  that 
out,  you  have  the  same  slope  out  to  what  corresponds  to  the  time  for  one  pass  of  a  satellite.  Then 
there  is  the  time  when  you  are  switching  to  another  satellite  and  you  get  a  different  type  of  behavior. 
I  agree  with  the  other  gentlemen  about  the  effect  of  the  periodicity,  but  this  is  not  that  effect. 


545/546 


